Nonlinear interference in a mean-field quantum model 
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Using similar nonlinear stationary mean- field models for Bose-Einstein Condensation of cold atoms 
and interacting electrons in a Quantum Dot, we propose to describe the original many-particle 
ground state as a one-particle statistical mixed state of the nonlinear eigenstates whose weights are 
provided by the eigenstate non-orthogonality. We search for physical grounds in the interpretation 
of our two main results, namely, quantum-classical nonlinear transition and interference between 
nonlinear eigenstates. 
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I. INTRODUCTION 

Commonly, quantum mechanical models of many in- 
teracting particles - linear models - are made computa- 
tionally tractable by transforming them into a mean-field 
single-particle models, at the cost of making the model 
nonlinear, but solvable [1]. Yet, solvable has a different 
meaning within two fields of condensed matter physics, 
the field of cold atoms in an external trap, and the field 
of electrons in a quantum dot. In the former field, the en- 
suing equation for the Bose-Einstein condensate (BEG), 
namely, the Gross-Pitaevskii equation (GPE), is often 
solved directly as a nonlinear differential equation [21 13] . 
In the latter case, for electrons confined in a quantum dot, 
the approximation applied for the Goulomb interaction 
between the electrons leads to the Hartree, the Hartree- 
Fock, or Kohn-Sham type of an equation. All these 
are nonlinear equations that are conventionally solved 
by an iteration scheme. These equations for the elec- 
trons have a similar nonlinear term as for the GPE, but 
the long range nature of the Goulomb interaction makes 
them more complex. Within each iteration, the mean- 
field potential is constructed with information about the 
charge distribution and/or the wave functions from for- 
mer iterations. The nonlinear Schrodinger-type equation 
is thus considered as a linear equation within each it- 
eration step and solved by means of methods from lin- 
ear algebra supplying a orthonormal set of eigenfunctions 
and corresponding eigenvalues. The nonlinear behavior 
enters the procedure once again when the eigenfunctions 
are used to calculate for the next step the new mean-field 
potential the electron is moving in. After the fulfillment 
of some criteria of convergence for the iterations, the end 
product is a solution in terms of orthonormal wave func- 
tions, presumably representing the single-electron states 
of the problem. 

In applications in quantum chemistry, it is often more 
appropriate to resort to functional basis that are not or- 
thogonal in order to reduce the size of the numerical ef- 
fort. A direct solution of the equations then leads to 
a general eigenvalue problem and solutions that are not 
orthonormal, but usually the orthonormality of the so- 



lutions is restored by a refined handling of the general 
eigenvalue problem [4 . 

For the fermion system (e.g. the electrons in a quan- 
tum dot), the hope is that the mean- field single-electron 
solutions that have been achieved by a self-consistent it- 
eration of the Schrodinger equation and the equation for 
the mean-field potential reflects in some reasonable ap- 
proximation properties of the original huge linear many- 
electron problem. Of course the individual single-electron 
states, wave functions, or orbit als are of limited value, 
but they can be used to construct more physically rel- 
evant entities like various response functions, the total 
energy, and the total charge distribution. Confidence in 
the method has come from comparison to numerical solu- 
tions of the corresponding truncated many-electron prob- 
lem for few electrons [5, 6j, and experimental results. 

For the boson system, namely, cold atoms in a trap, 
the ground state of the nonlinear GPE has been calcu- 
lated in order to gain information about properties of the 
BEG, and the GPE has been generalized to finite tem- 
peratures by the means of self-consistent Hartree- Fock or 
Bogolyubov approximations where special care has been 
taken to construct orthonormal states for the system [7] ; 
otherwise, reasons are given for neglecting the fact that 
a particular approximation does not preserve the orthog- 
onality [8 . 

After utilizing these methods for years in each of 
the mentioned subfields and comparing the different ap- 
proaches, we would like to draw attention of the reader to 
interesting open questions concerning the nonlinearity of 
the underlying equations. In particular, we show that the 
nonlinear mean-field single-particle ground state is not a 
pure state where all boson particles would condense, like 
in standard linear quantum theory. It is a mixed state 
that allows a small part (less than one percent) of the bo- 
son gas to populate higher nonlinear excited energy lev- 
els. We present a simple model to highlight our concerns. 
We are of course aware that these mean-field equations 
do only asymptotically describe the respective systems 
within an appropriate range of physical parameters like 
many other celebrated equations in physics, even though 
they have been found to reproduce the properties of the 
systems outside of the parameter range that the approx- 



imation itself can be justified for. 

In the summary section we reflect on the physical rel- 
evance and interpretation of our findings. 



II. THE NONLINEAR SCHRODINGER 
QUANTUM MODEL 

Prior to the discovery of the microscopic BCS theory 
of superconductivity as a Bose condensation of electrons 
paired by microscopic electron-lattice interactions ^, 
tentative links of superconductivity with the interaction 
between electrons and lattice vibrations had already been 
explored [10 [llj. Most interesting in the scope of the 
present work is Schafroth's early suggestion that charge- 
carrying bosons in a metal at low temperature actually 
constitute a gas of bound two-electron states [121 [E] 
which can then be described by a self-consistent nor- 
malized nonlinear Schrodinger-Poisson (SP) differential 
system [M]. At that time (1955), Schafroth just wanted 
to emphasize the role of long-range Coulomb interactions 
between bosons by use of mostly qualitative arguments. 
He assumed this charged boson gas to be defined by 
the stationary Schrodinger equation and proposed a one- 
dimensional (ID) SP model which considers the particles 
as moving in an effective potential V{x) related by Pois- 
son's equation to the local boson charge density with an 
additional uniform source term modeling both the charge 
density of the background (in order to make the whole 
system electrically neutral) and the non-condensed par- 
ticles. Evidently, this uniform source term in ID Pois- 
son equation can be regarded as the second derivative 
of an additional external harmonic confining potential. 
Consequently, this 1954 SP model was probably the first 
ab-initio mean-field nonlinear self-consistent attempt to 
describe the competition that occurs in a charged boson 
gas between external parabolic particle confinement and 
long-range internal particle-particle Coulomb repulsion. 

Even in ID, the numerical solution of the SP differen- 
tial system was not available at that time. Hence, the 
author considered the classical Thomas-Fermi (ThF) ap- 
proximation where both the effective potential V{x) and 
the boson gas wave function ^ were assumed uniform and 
respectively equal to the energy shift of the ground state 
due to Coulomb repulsion and to the wave function's av- 
eraged value L~^/^ over large regions of spatial extension 
L. These results led to the early conclusion that the oc- 
currence of superconductivity in metals is indeed due to 
the formation of some kind of resonant two-electron bo- 
son states at low temperature; a conjecture which soon 
became validated by the superconductivity BCS model as 
well as by Bogoliubov-type investigations of a Coulomb 
Bose gas [151 HI] • 

The present work revisits, by use of straightforward 
numerical routines, this pioneering SP differential model 
of a charged boson gas of particle mass M trapped 
in the radial axisymmetric parabolic potential V^(x) = 
V{r) = ^Mu'^r'^. As the simplest such example, let 



us mention opposite-spin Cooper-like electron pairs in 
quantum-dot helium [6 . Inasmuch as we wish to em- 
phasize interference-like quantum coherence effects due 
to nonlinearity, we extend our investigation to the neu- 
tral hard-sphere GPE boson gas. Indeed, the time- 
dependent interference between two BEC's has been in- 
vestigated in the one-dimensional case by exact solu- 
tions of the GPE obtained with methods of inverse scat- 
tering, thus displaying interference as an actual nonlin- 
ear interaction between GPE envelope solitons [17 . In 
this work, we also consider the complementary model 
of a single, harmonically-trapped, stationary BEG with 
discrete non-orthogonal (since the equation is nonlin- 
ear) interfering eigenstates |^i). Due to 2D radial ax- 
isymmetry, |^^) = l^ipi) ^ \m) is an eigenstate of the 
angular-momentum operator with the eigenvalue mh. As 
i'^il'^j) = {il)i\^l)j){mi\mj) ^ only if m^ = rrij = 0, 
it is sufficient to consider the zero-angular-momentum 
states 777, = in order to display peculiar physical prop- 
erties related to Bose eigenstate non-orthogonality (in 
contrast with, e.g., BEC vortex- nucleation where the 
777, = 1 nonlinear eigenstate, actually orthogonal to the 
777, = ground state, plays the major role [18 ). The 
discrete radial one-particle Bose eigenfunction '^^(r) is 
defined in the mean-field approximation by the station- 
ary Schrodinger equation related to the corresponding 
chemical-potential eigenvalue jUi. In scaled form Ui{r) ex 
^l^i{r), it reads: 
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where the dot stands for derivation with respect to the 
radius X meas ured in units of the characteristic length 
Iq = y/h/2Muj and the tilde superscript labels energy in 
units of hw. The dimensionless particle-particle interac- 
tion energy ^i(r) per particle introduces into the system 
nonlinearity which is appropriately scaled to unity. In- 
deed, defining 



^iW = V i^^^W' 



(2) 



where the dimensionless parameter A/i is self-consistently 
given by the solution of Eq. ^ according to 



/ u'^iXdX=JVi, 
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(3) 



we obtain the necessary 2D axisymmetric normalization 
condition 



/ \<i/ifd^x= f ^p^2TTrdr = l, 



(4) 



together with the scahng to unity of the nonhnear coef- 
ficient for both the GPE: 



*i = ul 



(5) 



and the SP differential system 
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Equation ^ defines the contact interactions of the hard- 
sphere BEC while the 2D Poisson Eq. (|6| defines the 
long-range Coulomb interactions of the charged Bose gas 
(e.g. quantum-dot helium 16 ). Practically, the value of 
the dimensionless norm ([3| is given by the experimen- 
tal conditions. We have Afi = 4.N{as/Lz) for the BEC 
defined by its scattering length a^ and extension Lz in 
the axial z-direction that contains N particles in the zth 
nonlinear eigenstate [18 . Specifically, A/i = 30.29 for 
the ground state in Paris-ENS' large laser-beam stirrer 
experiment pHST] while Mq = 187.35 in Boulder- JILA's 
rotating normal-cloud experiment [22, 23 . These high 
values JVi :^ 1 of the nonlinear parameter are due to 
the number of atoms in the trap which is quite large 
(typically of order 10^~^). They allow the ground state 
to be approximated by its localized "negative-curvature" 
parabolic Thomas- Fermi profile [2j |3] and define in the 
present work the so-called "classical- nonlinear" regime. 
On the other hand, the experimental parameters defining 
quantum-dot helium yields much lower values A/i ~ 2 — 3 
for the ground state nonlinearity [24 ; they would corre- 
spond in our study to the "quantum-nonlinear" regime 
where interference effects occur. 



III. DISCUSSION 

The interest in nonlinear quantum eigenstates can be 
illustrated for both the charged Bose gas defined by Eq. 
(|6| and the GPE system defined by Eq. (5|. In the 
former case, the nonlinear eigenstates can be regarded 
as unperturbed although they take into account, in ad- 
dition to the particles' external parabolic confinement, 
the (usually quite important) long-range Coulomb in- 
teractions ^. In the later case, nonlinear quantum 
eigenstates in the physical description of very-many- 
particle stationary BEC's are unavoidable. Indeed, there 
is practically no other choice than solving the nonlin- 
ear GPE (Tl) and ([s]). Therefore we numerically con- 
sider the discrete real-valued radial-symmetric normal- 
ized eigenstates Ui of Eqs. (fTl) and ([5|-(p| which are 
by choice non-orthogonal since we restrict ourselves to 
rrii = zero-angular- momentum s states [26 . Specifi- 
cally, we consider the two first such eigenstates \u^) and 
\u^) whose superscript labels refer to their number of 
huu eigenquanta in the limit of vanishing nonlinearity, i.e. 
for the 2D linear parabolic system [6]. Their normalized 
inner product {u^\u'^) ^ yields the statistical weight 
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\{u^\u^)\^ = {u^\u^)^ + that defines 
the mixed ground state whose appropriate description in 
terms of the nonlinear density matrix p is 
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As a matter of fact, the standard use of a density matrix 
demands an orthonormal basis of eigenstates. Therefore, 
Eq. ^\ should be understood as an extension of this con- 
cept to the nonlinear case where the eigenstates are not 
orthogonal. This generalization is all the more natural 
as the deviation from orthogonality is small (see below 
Fig. fl]) and it yields a transparent physical interpreta- 
tion. Indeed, Eq. ^\ suggests that the probability per 
particle for the Bose gas to be in the nonlinear excited 
state \v?) is but its probability 1/(1 + ur"'^) to be in the 
nonlinear ground state \u^) multiplied by the mere tran- 
sition probability w^^ from \u^) to \i^). However this 
last expression must be taken in a somewhat loose sense 
since the present time-independent description is, strictly 
speaking, incompatible with the concept of a quantum 
transition (we will return to this discussion further be- 
low). 

The paradoxical property, compared with standard lin- 
ear theory, is that as the system is in the mixed state 
described by Eq. ([t]), it nevertheless allows interference 
between the two eigenstates \u^) and \i^) to occur, ac- 
cording to the following theorem: 
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The subscripts in ^\ define the Hilbertian matrix el- 
ements of the particle-particle interaction potential $* 
corresponding to \u^) (see Eq. ([I])). These elements have 
been calculated by use of the nonlinear eigenfunctions u^ 
while Pi are those nonlinear eigenvalues (or equivalently 
chemical potentials) which respectively define the ix*'s. 
Equation (Is]) is a direct consequence of the Hermiticity 
of the Laplacian operator in the Schrodinger equation 
^ [27 . The first term on the r.h.s. of Eq. ^ defines 
the probability amplitude for the system in the nonlinear 
eigenstate \u^) to be also in the nonlinear eigenstate |ii^), 
due to the interaction potential <I>^ defined by the proba- 
bility density \u^\^ = {u^)^ through Eqs. ^ or ([6|, while 
the second term defines the reverse process, namely, the 
probability amplitude for the system in \u^) to be also 
in \v}) as a consequence of the interaction potential $^ 
defined by [u^Y . Equation ([8| is exact. Therefore no 
perturbative-like ordering in <I>* is needed although, of 
course, a straightforward time-independent perturbation 
scheme that considers <l>* as perturbation in its respec- 
tive Schrodinger equation ([l]) allows one to recover it (for 
instance when nonlinearity is weak). This point will be 
further developed in the next part. 

The two amplitudes in Eq. ([s]) interfere because 
the corresponding 1^3 processes are indistinguish- 
able in the build-up of the probability amplitude 
l\Vi^\u^). Therefore this amplitude actually defines the 
nonlinear quantum coherence in our two-state system. 
The interference pattern is increased by the a; ^ pro- 
gressive flattening of the parabolic trap in the case of 
the confined charged Bose gas and by the increase of the 
particle number in the GPE system; or, equivalently for 



both of them, by a progressive increase of the nonhnear- 
ity parameter M in the system, as defined by Eq. (|3| . 

Figure [l] illustrates in the case of 2D axisymmetric 
quantum dot helium the remarkable behavior of the nor- 
malized mixed-state statistical weight 
\{u^\u^)\^ = {u'\u^)^ where 
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The two eigenstates u^^"^ must obviously correspond to 



In the case of long-range Coulomb particle-particle in- 
teractions described by Eq. (|6|, there exists a maximum 
amplitude uq for both nonlinear eigenstates u^ and u^ 
that is clearly visible on Fig. 2 when JV varies from 10~^ 
{uq ~ 0.1) to 10^ {uo ~ 1). For such high values of A/", 
one reaches the classical asymptotic ThF regime defined 
by neglecting the quantum kinetic Laplacian derivative 
terms in Eq. ^ [14 . Then <I^(X) - /i^ - X^/4 and Eq. 
([6| yield the common limit u^{X) = u^{X) = utf = 1- 
The respective initial conditions of the two modes spi- 
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FIG. 1: The quantum-dot helium interference pattern defined 
by the square probability amplitude w^^ = (u^\u^)^, as com- 
pared with its sin^ approximation defined by Eq. (Il0| (thin 
continuous line). There occurs a clearly visible threshold at 
A/" ~ 6 where the system bifurcates from its interfering quan- 
tum regime towards its Thomas-Fermi classical one. 

the same physical system. For the GPE defined by Eq. 
(pi, this means that their particle number should be the 
same and hence A/" = A/i = A/3. The case of the con- 
fined charged Bose gas defined by Eq. (|6| demands an 
additional condition which states that the external trap 
parabolicity uj is identical for both eigenstates u'^. Never- 
theless it also yields A/" = A/i ~ A/3 [24] . We numerically 
obtain the simple interference pattern 
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for A/" < 6 and 7 = 2.29 10"^ It displays a re- 
markable nonlinear resonance where the quantum co- 
herence maximum is reached about the particular value 
A/" = A/"* ~ 3. For quantum-dot helium, the only free 
physical parameter is the trap parabolicity uj. Therefore 
we have A/" = A/'(a;) and the resonance occurs for the 
particular parabolic trap profile V{r) = \Muj*r'^ where 
cj* = a;(A/'*). Numerical simulations in this 2D axisym- 
metric model show that ficj* ~ 0.14 e where e = Me^ jT? 
is the effective atomic energy unit [24 . Hence ficj* ~ 3.80 
eV for electrons in vacuum since e = 27.21 eV while 
ficj* ~ 1.66 meV in the case of GaAs quantum-dot helium 
where e = 11.86 meV. 




FIG. 2: The convergence in the {Co = /i - l>(0)} vs {1x0 = 
1^(0)} initial-condition phase space of the two- level SP non- 
linear system defined by ix^ (lower red curve) and i^ (up- 
per blue curve) towards the common fixed point u^ — \ and 
Co = when the nonlinearity parameter increases from the 
linear regime M ~ 10~^ to the classical Thomas-Fermi one 
defined hy M ^ 10^. The corresponding common asymptotic 
eigenstate profile utf (X^ for these two first ?7i = nonlinear 
eigenmodes is the uniform one utf = 1- 

ral down towards the same Thomas- Fermi fixed point for 
A/" ^ 00. At the same time, the corresponding eigenstate 
profiles ui and us start increasing their width instead of 
their amplitude: see Fig. [3J There is a progressive merg- 
ing of ui and us into the single asymptotic uniform mode 
Utf defined by the fixed point {1, 0} in Fig.|2] Quantum 
discreteness disappears for A/" ^ 1, which is the hallmark 
of the classical regime. 

Quite different is the convergence of the GPE sys- 
tem ([5]) towards the ThF classical regime still defined 
by l>^(X) - /ii - XV4, as shown by Fig. ^ Then 
we indeed have lim^^^^oo C'o = lim7vr^oo[M ~ ^(0)] = 
lim^v'^ooi/^ ~ '^0] — O5 t)ut the nonlinear eigenstate am- 
plitudes Uq grow without limit while thei r correspo nd- 
ing profiles converge towards utf{X) = ^s/jl — X2/4, as 
shown by Eqs. [T] and [5] when the Laplacian derivative 
terms are neglected. There is no bifurcation from the 
amplitude-growing regime to the width- growing one, like 
in the SP case displayed by Figs. [2] and |3] It can be said 
that, due to GPE particle-particle contact interactions 
defined by Eq. ([5|, the bosons pile up in the trap rather 
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FIG. 3: The nonlinear SP eigenstate profiles u^ (continuous) 
and u^ (dashed) for increasing values A/" < 8 of the nonlinear- 
ity. Note the last u^ and u^ profiles where the width — instead 
of the amplitude — starts increasing. 



FIG. 5: The interference pattern defined by the square prob- 
ability amplitude w^^ = {u'^\u^)^ in the case of the GPE 
system. The threshold for the quantum-classical transition 
occurs at J\f ^ 10. 




FIG. 4: The convergence in the {Co = /i - $(0) = jl - ul} 
vs {uo = u(0)} initial-condition phase space of the three-level 
GPE system defined by its ?7i = nonlinear states u^ (lower 
blue curve), u^ (middle pink curve) and u^ (upper red curve) 
towards the common fixed point uq = oo and Co = when 
the nonlinearity parameter increases from the linear regime 
JV ^ 10~^ towards the classical Thomas-Fermi one JV ^ 
1. The corresponding common asymptotic eigenstate profile 
utf{X) for these nonlinear eigenmodes is defined by utf = 
y^jl — X'^/4 as shown by Eqs. Ill and l5| when the Laplacian 
derivative terms are neglected. 



than spread out. The corresponding nonlinear resonance 
defined by Eq. ^ with A/" = A/i = A/3 looks also quite 
differently, as shown by Fig. [5] 



IV. CONCLUSION 

The nonlinear eigenstates \i) that have been defined 
and discussed in the present work are time-independent. 



Since the nonlinear term in Eq. (IT]) is a (square) mod- 
ulus in the corresponding time-dependent Schrodinger 
equation, its eigenstates ^* ex u'^ e*^^^/^ do still have the 
standard time dependence related to the eigenvalue /i^ 
of \i). A conceptual problem arises when one considers 
the interaction between such states. Indeed, due to their 
non-orthogonality defined by Eq. ([9| and illustrated by 
Figs. [1] and [5| one cannot state any more that the sys- 
tem is in the pure stationary eigenstate \i). It has the 
probability amplitude {j\i) defined by Eq. ([8| to be also 
in some other nonlinear state \j). The statistical weight 
w'^^ = {j\i)'^ then defines the corresponding mixed state 
in accordance with Eq. ^. Actually w^^ is quite small: 
typically less than 1 % (cf. Eq. (10)). This might explain 
the quite acceptable results obtained in GPE systems 
when one neglects w'^^ and assumes that all the boson 
particles have condensed in the ground state. Neverthe- 
less, this statement is formally wrong and might yield 
a weak measurable departure from the expected GPE 
ground state particle density in accurate BEG experi- 
ments. 

It is tempting to relate w'^^ to a mere transition prob- 
ability. Again, the result is unexpected. Indeed, assume 
that the nonlinearity Af = Afi = Afj is small. There- 
fore each interaction potential <l>*'-^ may be regarded as 
a perturbation in its corresponding Schrodinger equation 
([1]). Using standard time-independent perturbation the- 
ory, one obtains ^lj/{Ei—Ej) (resp. ^^j-/{Ej—Ei)) as the 
probability amplitude for the system being in ith (resp. 
jth) linear eigenstate of energy Ei (resp. Ej) to be also 
in the jth (resp. ith) linear eigenstate of energy Ej (resp. 
Ei). One recognizes in the two terms of Eq. ([8| the sim- 
ple extrapolation of these probability amplitudes to any 
value of the perturbation (i.e. of the nonlinearity) while 
the linear energy eigenvalues Eij become the correspond- 
ing (nonlinear) chemical-potential fiij. However, due to 



nonlinear quantum coherence, there is no simple rela- 
tionship between each such probability amplitude and a 
possible transition probability. Indeed, for the {|z), \j)} 
two-state system, exact diagonalization of the perturbed 
Hamiltonian in the time-dependent Schrodinger equa- 
tion is possible. It yields in first-order with respect to 
^ the so-called Fermi golden rule transition probability 

\i) ^ li) 
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or reverse 
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Since Eq. (11) is exact, one can average 



it for t > h/{Ei - Ej) and obtain Vi 



^j 
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Ej) in the lowest order in <l>^; which, apart from the fac- 
tor 2, is indeed the square time-independent probability 
amplitude that has been previously obtained. However, 

{j\iy is not 



i(^wi+^i- 



due to Eq. rt8|), the statistical weight w'^^ 

merely equal to the extrapolation ofP = 

to any value of the nonlinearity parameter J\f. There 

remains twice the probability amplitude product which 

lowers P by more than one order of magnitude. 

By focusing our attention on the nonlinear charac- 
ter of two simple cases, the present paper actually aims 
to be introductory in the field concerning the reduction 
of the original 3N-dimensional ground-state wave vector 
defining a N-particle linear quantum system to the 3D 
one-particle mean-field nonlinear Schrodinger description 
based on the mixed ground state given by Eq. ^. The 
conclusion that the nonlinear Schrodinger equation, like 
many other similar ones, is only intended as a recipe 
to obtain a solution within an iteration scheme where 
the problem is considered linear at each iteration step 
lies close at hand. However, it remains an open ques- 



tion whether the direct solution of the nonlinear problem 
yields properties that are not obtainable in the iteration 
scheme, and indeed, if such inherent nonlinear properties 
have a physical relevance. While the (mostly numerical) 
investigations of the time-dependent GPE mainly con- 
cern linear surface excitations of the condensate about 
its stationary asymptotic ThF ground state profile, true 
time-dependent nonlinear structures such as the breath- 
ing monopole oscillation mode have been pointed out [2] . 
These structures are actually solitons, like in the case 
of two-BEC time-dependent interferences [17]. We re- 
call that solitons provide a paradigm example of a true 
nonlinear phenomenon that can never be described as 
a convergent limit of a linear iterative process. Such 
seems also the case for the quantum properties related 
to non-orthogonality between the one-particle mean-field 
nonlinear eigenstates. Therefore any attempt to investi- 
gate the time-dependent relationship between linear and 
nonlinear eigenstates of a given quantum system — like 
the recent attempt to both display particular nonlinear 
eigenfunctions as global attractors for all finite-energy so- 
lutions and describe quantum transitions between them 
[28] — are welcome. However, only a comparison with ex- 
periments and the solutions of the original linear many- 
particle problem can shed a definite light on these ques- 
tions. 
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